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Abstract 

The definition of a mobilized system and its efficiency are introduced. The existence 
of an optimal (maximally efficient) system is proved by an application of Young measures 
and compensated compactness. 



1 Introduction 

1.1 Motivation 



Which one is the exception- a motor boat, a car, an airplane or a submarine? The immediate 
(perhaps, after a short contemplation) answer is, of course, the car. Indeed, a car is the only 
vehicle in the above list which "defies" the law of conservation of linear momentum. The 
others create an opposite stream in the medium in which they move (air, water) which, by 
conservation of linear momentum, must be compensated by the motion of the vehicle. Take 
away the friction of these vehicles and the medium, and they will perform better (faster, more 
efficient). 

The principle behind the motion of a car is different. It is moving (that is, shifting from 
rest to cruise velocity) because of the friction created by the contact of its tires with the road. 
Take away the friction and the car will not be able to move at all, independently of how 
hard you push the gas pedal. You will not be able to stop either, if your friend, unwisely 
attempting to help, gave you an initial push on a frictionless road. 

So, what is a car? If we strip it off the non-essential components (radio, GPS, the fashioned 
seat covers etc.), it is a collection of components which can move with respect to each other 
under the preassign constraints caused by the mechanical structure. Unless you change gear 
or the pressure you apply on the gas pedal, the motion of these inner part can be assumed 
to be periodic (again, with respect to the frame of reference of the car itself). This concept 
has a lot in common with the subject of optimal locomotion of a swimmer in a Stokes flow 
- where the motion is due to a periodic change of shape of the swimmer in the absence of 
inertia. The description below is, in a sense, the mechanical analog of the swimmer model, 
see [AGK] and the references listed there. Another series of publications which seems to 
be related to the present discussion concerns molecular motors and the flashing rachet. See 
[CHK], [CKK], DKK] and ref. therein. 

Note that in the case of micro-swimmers and molecular motors, it seems that there is 
no intuitive way to predict the direction and velocity of the swimmer (res. motor) from its 
periodic motion. A car, at a first glance, is different. However, the abstraction of a car which 
we consider in this paper (and call a "mechanical system") is, in a way, a generalization of 
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the concept which contains cars, microswimers, molecular motors and, perhaps, many other 
objects whose dynamics are not inertial in nature. 

1.2 Objectives and outline of results 

A mechanical system is represented by a Lagrangian L = L{x, t). It is T— periodic in the time 
t, and x{t) stands for the position of a reference point (say, the center of mass of the system). 
Such a system is called "mobilized" if the global minimizer x{t) of the action L(x,t)dt 
is not periodic, i.e x{T) ^ a;(0). In section |21 we attempt to justify this model and explain 
the reason why the global minimum of the action represents the asymptotic motion of the 
system under friction. In addition, we introduce a reasonable definition of efficiency, in terms 
of a relation involving the speed v := |x(T) — x(0)|/r due to the action minimizer, and the 
minimal action D = L{x,t)dt itself. The number D stands for some indication of the 
energy dissipated (or invested) per period, so a more efficient system means larger v and 
smaller D. We scale the efficiency function of a mechanical system L in such a way that 
< < 1. Then, we ask if either there is a (theoretical) possibility to achieve the "most 
efficient" system L whose efficiency e := e-^ < 1, or there is a sequence of systems L„ whose 
efficiencies eL„ — > e < 1, but the "ideal", most efficient system L does not exist. 

In section IHl we consider a special case of mechanical systems, where L is a homogeneous 
function of x, and prove the ffist alternative: There is a "best" (most efficient) mechanical 
system and its efficiency e is always smaller than 1. 

2 Description of the model 

Let us attempt to build a mathematical caricature the " car" , composed of a finite number of 
"parts". To wit, assume it is composed of a collection of N points of respective masses mj 
executing orbits Xi = Xi{t) on the line R, z = 1, . . . n, so that 



where v is the effective velocity of the car (with respect to a reference frame attached to the 
road) and T is the period of one cycle. 

We assume mj = 1. The orbit Xi{t) is given with respect to a fixed frame of reference 
(attached to the road). We may write it as 



where yi is the orbit of the corresponding point with respect to a reference attached to the 
car and x = x{t) is the orbit of the whole with respect to a reference attached to the 

road. In this representation yj is a periodic orbit, representing the motion of a part forced 
by, say, the internal combustion of the engine. 

The simplest model for the motion of x{t) is the linear forced system 



Xi{t + T) =Xi{t)+Tv 




Xi{t) = yi{t) + x{t) 



x + [3x = F{t) 
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where /? > is the friction coefficient and F is the total forcing acted on the car by the 
motion of its inner parts. It is the sum of the forces F := /j, where fi, the force appHed 
by the i part, is defined in terms of yj as: 



fi{t) = -mi{yi + (3yi) 



Now, define 



N 



(2.3) 



(2.4) 



i=l 



Then ()2.21 12. 3 j) can be summarized as the Euler-Lagrange equation corresponding to the 
Lagrangian 



£W(x,t) :=e^*/:^^>(x,t) 



where 



L(s) := |s|V2 . 



(2.5) 
(2.6) 



Indeed, the Euler-Lagrange equation associated with C^^^ under (|2.()j) is 



= e' 



N 



N 



^ mi{x + 2/i) + /3 ^ mi{x + iji 



which implies (|2.2|) via ()2.3|) and the condition ^rrii = 1. 

The energy dissipated per cycle for an orbit yi, . . . y^r, x is given by (3D where 



TV 



mi 



i=l 



dt= c^y^x,t)dt . 



(2.7) 



We now generalize (|2.4I2.7() into: 



Definition 2.1. A mechanical system is determined by a forced orbit composed of N 

periodic functions y = {yi{t), . . . yN{t)} in terms of the Lagrangian as given in \2. 5|] 
and a convex function L generalizing h2.b]) . 

The moment associated with C^^^ is denoted hy p = cS^^ . The Euler-Lagrange equation 
associated with (|2.5|) is 

p + = ^ p{t) 

so the asymptotic motion of a mechanical system is determined by the orbit p{t) = 0. 
To elaborate, let 



n^y}{p,t) = sup 

c 



p.C-/:{y}{C,t) 



(2.8) 



be the Hamiltonian associated with the Lagrangian £.^y\ The equation of motion corre- 
sponding to the non-dissipative dynamics is given by 



x(t) =7t:M(A,t) 
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where A stands for the constant momentum p. If a friction is apphed, then x = x{t) is 
determined by p = 0, namely 

x{t) = n}y^o,t) 

is the asymptotic motion of the system. It is a global minimizer of the action determined by 
the Lagrangian £^y\ Note that 

T 



min 

x=x{t) Jq Jq 



Let 

»(0) = ijf'H<.>(0,tMt.f(azfi(£), 

where x{t) is the global minimizer of the action. 

Definition 2.2. A mechanical system for which v{0) ^ is called a mobilized system. 
The first result is somewhat disappointing: 

Theorem 1. If L is a quadratic function \2.b\) . then the system is not mobilized . 
Proof. 

N N 

c^y^x, t) = -Y,m \m + =^ c^f = E (^*(*) + • 

i=\ 1=1 

In particular, = implies x = — X^i^^w-jyi- Since yi are periodic by definition, then 
/(f xdt = v{0) = 0. □ 

We now generalize the energy dissipation in the linear case (|'2.7j) for general mechanical 
systems. We shall denote by D the minimal action, and refer to it as the energy dissipated 
along a cycle: 



D = min 

x=x{t) 



^j^ ciy}{i,t)dt = -^l^ n^yHo,t)dt. 



Next, we define the efficiency of a mobilized system y := yi, . . .y^- This should indicate 
the ratio of the distance transversed per cycle to the dissipated energy. The right scaling for 
this turns out to be 



eL{yi,---yN) ■■= — = — = — ^ ^ rT^,s X. — ■ f^^-^^) 
D T-i /q n^y}{o,t)dt 

In fact, it can be proven that 

Lemma 2.1. If L is a convex function, then for any mechanical system composed of N 
periodic orbits yi{t) = yi{T -\- 1) , i = 1, . . . N , 

< eL{yi,...yN) < 1 • 
//, moreover, L is strictly convex, then cl < 1. 
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Proof. Let x(t) be the any orbit. Then, 

rT N 



T-^ [ C^y\±,t)dt = T-^y^mi f L{y,{t) + ±{t))dt . 
Jo Jo 

By Jensen's inequahty, the normahzation condition ^rrii = 1, the periodicity of yt and the 
convexity of L: 



T'^Y^rmJ^ L{yi{t) + ±{t))>L(T-^J^ J2miiyi + ±)dt^ =l(^T-^ ±d?j , (2.11 



so < < 1. Now, if L is strictly convex then the equahty in (|2.1ip holds if and only if 
N = 1 or yi = yj for all 1 < i, j < N . But, in the later cases, the optimal orbit x is evidently 
equal to any component yj, so it is a periodic function and v{0) = 0. Hence the system is not 
mobilized and efficiency is not defined (or cl = altogether). □ 

Assume now that a convex Lagrangian L is given, as well as A*" G N and {mi, . . .tun} G 
Ei^i = 1. Let 

AL:={2/ = (yi,...yiv)GCi([0,l];M^) ;y(0) = y(l) .} (2.12) 

Let 

Ci := supeL(y) . (2.13) 
yeA 

We know that bl < 1. The intriguing questions are 

i) Is CL = 1? 

ii) licL < 1, can the supremum in ()2.13() be achieved in some sense? 

We try to answer these questions in the special case of homogeneous Lagrangians. 



3 Homogeneous Lagrangians 

To fix the idea, let us concentrate on the case L(s) = \s\" for some a > 1. We shall further 
assume a unit period T = 1. 

Let ^ 

H" := |x = x{t) : [0, 1] ^ M , < t < 1 ; ^ < ooj . (3.1) 

Let also 

:= {y e {H'^f ,y{0) = y{l)] . (3.2) 
Given rrii = m £ (0, 1), rrii = I and y = (yi, . . . yN) £ Aq- consider the Lagrangian 

N 

C{x,y) ■.= '^mi\yi + xl'' . (3.3) 
1=1 
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Define also for any such y and p £ M^, 

n{p,y):= snv {p-i-C{i,y)} , (3.4) 



and 



Also, define: 

u{y) : 

and 



H^(p) := [' n{p,yit))dt (3.5) 
Jo 

dn{p,y) 



dp 



(3.6) 



p=0 



r-l 

< u >y.= I u{y{t))dt . (3.7) 







We obtain via (l;-!.4l;-i.7ll . (tTTUl) and (tTTTl) 

ea{y) ■= ; = sup e^(?/) . (3.8) 

One may wonder if, under the homogeneity condition, the mechanical system is mobilized, 
i.e > 0. The first result we claim is that this condition is equivalent to a property of the 
function u as defined in (|3.6jl . namely 

Lemma 3.1. The homogeneous system is mobilized if and only if u is not a linear function. 

Proof. If u is a linear function, then for any y £ A^, J u{y)dt = by definition of A^. 
Indeed, the condition 7/(0) = y{l) is equivalent to y{t)dt = 0. 

Conversely, if u is not a linear function, then there exists € Ao- for which u (y) dt ^ 0. 
This implies e^iy) > and, in particular, the system is mobilized. □ 

The condition of non-linearity of u is rather delicate. In fact. 
Lemma 3.2. // either a = 2 or N = 2 then u is a linear function. 

Proof. We already know that a = 2 (quadratic Lagrangian) is not mobilized for any N £ N 
by Theorem^ 

By the homogeneity of L and the definition of u we observe that u satisfies the pair of 
symmetries: 

VCeM^ , e = (!,...!), and AG M , u(C + Ae) = n(C) + A 



Va G M+ , u{aC) = au{C) (3.9) 

We conclude, therefore, that u{yi,y2) = f{yi — ^2) + Vi holds for some function / of a single 
variable. From the second equality of (|3.9|) it follows that 

f{a{yi - 2/2)) + ayi = af{yi - 7/2) + ayi =^ af{C) = /("O 
for any a G M and ^ G M. Hence, / is linear and so is k. □ 
The main result of this paper is: 



6 



3.1 Main result 



Theorem 2. // the homogeneous mechanical system \3. 5)) is mobilized then there exists a 
maximizer of in the set Ao- , and Co- < 1 • 

Remark: Lemma 13.21 implies that N > 2 and a ^ 2 are necessary for the condition of 
Theorem [5J We conjecture that these conditions are also sufficient. In any case, it is not 
difficult to construct examples for mobilized homogeneous systems. For example, take a = 3, 
N = 3 and mi = m2 = = 1/3. The function u = u{yi,y2-,yz) can be readily calculated as 
the root of a quadratic equation whose coefficients are linear functions of yi. The discriminant, 
however, is not a complete square, so u is not linear. 

In the rest of this section we prove Theorem |2I 

From ()3.9|) we obtain that "^(0, y) is a— homogeneous, that is 

a G M , 7^(0, ay) = |a|'^W(0, y) 

as well as 

< U >ay= Oi < U >y . 

In particular 

ea{ay) = e^{y) , Va G M, G . (3.10) 

In addition, eij{y) is clearly invariant under diagonal shifts y y{t) + I3{t)e where e = 
(1, . . . 1) G M^. Define now 

A° := {y = {yi, . . . yjv) G ; yi = .} 

and 

5. = |yGA0; l\y{t)Wdt = l^ ; = !^y G ; |(y(t))r < l| . (3.11) 

By the scaling (|3.inj) and the diagonal shift invariance we conclude that 

Ca := sup Caiy) = sup e^iy) = sup e^iy) ■ (3.12) 

yeAa y&Sa yeBa-{o} 

Let now yj be a maximizing sequence of e^ in S^. There is a weak limit y^ G B„ of this 
sequence. The inequality 

hm (0) < H^- (0) 

holds since H^(0) is upper-semi-continuous, but we do not have the same claim for < u >y. 
So, we cannot prove that vecy oo is a maximizer of e^. 

Another problem is that we may have y^ = 0, so eo-d/oo) is not defined at all. As 
an example, let y G A^ and assume that < u >y^ 0. This, in particular, implies that 
x{t) := u{y) = x{t) + \t where x is a periodic function and A 7^ 0. If we replace y 
by yj = Ujif) '■= j^^yU^) for j G N, using the periodicity of y to define yj on (0, 1), the 
following claims are straightforward: 
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a) yj G Ao- for any j G N. 

b) linij^^oo yj = weakly. 

c) Xj = r^x{jt) = j'^x{jt) + \t = J^U (jjj). 

d) < u >y. = A for any j € N. 

In particular, we find out that e^iyj) = eo-(y), while is not defined for the weak limit 
limj^oo yj = 0. 

3.2 Relaxation 

To overcome the last difficulty we shall extend the definition ()3.5II3.7() as follows: Let 

:= {Probability Borel measures on M^} 

Given v eV^ , let 

H^(p) := f n{p,v)v{dv) (3.13) 

and 

< u >u:= I u{v)v{dv) . (3.14) 
The generalization of (|3.8|) is given by 



|< U >y\ _ \< U>1 



We shall further extend the definitions of Aq- and as follows: 

■= \ v = v{dv) eV^ ; / |t;|'^i^(dt>) < cx) ; / Viu{dv) = ,1 < i < n\ (3.16) 

a|| := {z^ G Ao- ; = 6v-^fi{dv2, . . . dv^) where ji G (3.17) 

Lemma 3.3. For any y G Aq- (^res. y G A[^J i/iere exists v ^ tv^ (res. v G A^^j so t/iai 
ea{v) = efj{y). Conversely, for any v d A.^ (res. v G A^J i/iere exists y G Ag- (^res. y G A[^J 
so t/iai eo-(z^) = eaiy). 

Remark: The measure v associated with y is related to Young measure. In general, however, 
Young measures are used to study the oscillatory behavior of a weak limit of L°° sequences 
(see, e.g., [E]). 

Proof. For the first part, define ^{dv) = JqYIi ^Vi-yi{t)- To elaborate, the measure u corre- 
sponding to y is obtained by its application on test functions cp £ Cq 

(f){v)u{dv) = [ (p{y{t))dt . 
Jo 
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If y G Act then the above equahty also extend to <^(y) = y and (t){y) = \y\^ ■ In particular 

[ vu{dv) = [ y{t)dt = ; / \v\''v{dv) = [ |yr it)dt < oo . 
Jr^ Jo jr^ Jo 

Finally, the equalities 

H^(0) = H^(0) ; < u >y=< u >y (3.18) 

hold under this identification. 

For the second part we use Theorem 2.1 of [Am] to observe the following: 
For any such v there exists a Borel function T : [0, 1] — > which push forward the Lebesgue 
measure dt on [0, 1] to v. That is, for any test function cj) £ Cq 



f (t){T{t))dt= [ ^{v)v{dv) 
JR^ 



In fact, Theorem 2.1 of [Am] claims the equality between the infimum of Monge and the 
minimum of the Kantorowich transport plan of probability measures to /U2 on M^, provided 
/ii has no atoms. This implies, in particular, that the set of Borel mappings transporting 
^1 to /i2 is not empty. In our case we use this result under the identification of with 
the Lebesgue measure on [0, 1] (considered as a Hausdorff measure in for the embedded 
interval) and ^2 with v. 

So, set y{t) := f^T{s)ds. Is is absolutely continues function satisfying y = T a.e. The 
definition (|3.16|) implies also that Jq \y\^ = J \v\'^i'{dv) < 00, as well as /q ydt = J vv{dv) = 
0, which yields the periodicity of y. The equality (|3.18|) holds under this identification as 
well. □ 



3.3 Proof of Theorem [H 

Define now 

TV 

i){v) := \u{v)\ + '^\vi + u{v)\ , V = {{),V2,...vn) ■ (3.19) 

i=2 

We shall summertime some properties of ij) which will be needed later: 
Lemma 3.4. There exists A> Q so that 

A~^\v\<\%l){v)\<A\v\ "iv = {{),V2,...vn) . (3.20) 
In addition, for any € A^, 

A-^ [ \'ilj{v)\''u{dv) < -H^(0) < A [ \iP{v)\''u{dv). (3.21) 
Jrn JrJv 

Proof. The estimate 1)3. 20() follows from the homogeneity of u, namely u{av) = au{v), as 
well as from the evident property n(0, v, . . .v) =0 <^=^ v = 0. The estimate 1)3. 21[) follows 
from (ESI 031) and ffTH^ . as weh as (jT^ . □ 
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Given 1 < q < a, let be the unit sphere (res. the unit ball) defined by 
Sl = LeAl; [ \v\'^u{dv) = l\ ; X = L gA^ ; [ \v\'^u{dv) < 1 } {3.22) 



The analogous of (|H.12j) also holds due to the scaling invariance 
Lemma 3.5. Cg- = sup^^^g eo-(i^) • 

Proof. We only have to show that any G A^- (res. u £ A^) can be transformed into 
eSa (res. i> G 5°) so that e^(zy) = e^(z>). Let u{dv) = l3^u{pdv) for /3 > 0. The 
homogeneity properties (|3.1U|) imply that, indeed, eo-(i>) = eo-(^') for any such 13. In addition, 
j vu{dv) = Jvi/{dv) = if e Ao-. However, / |v|^z>((it>) = P^'^ J \v\'^i^{dv), so i> e 5o- if 

p= {Jlvl^uidv))^^'. □ 

Let Vj be a maximizing sequence of Cq- in 5"!^. Since the q > 1 moments of Uj are uniformly 
bounded, this sequence is compact (tight) in the weak topology of measures. Let i/qo be the 
weak limit of z^j. Since \u{v)\ < A\v\ for some ^ > it follows that the sequence of (signed) 
measures u[v)uj[dv) is tight as well, and that 

lim u{v)vj{dv) = u{v)voo{dv) . (3.23) 

On the other hand, it is not a-priori evident that Voo G S^a^ since the sequence \v\'^Vj is not 
necessarily tight, so we only know 

\v\'ivoo{dv) < 1 . (3.24) 



We claim, however, that, in fact, the a moments of Uj are uniformly bounded. For, if 
/ \v\'^Vj{dv) oo, then by 1)3.201 IS^^T]) . also — H'^j (O) — > oo. Since we also now that < u 
are uniformly bounded, then 

= -^m " ' ' 

contradicting the assumption that fj is a maximum sequence for e^- 

Since q < a hy assumption it follows that \v\^i'j is a tight sequence as well, so there is, 
in fact, an equality in ()3.24)) . In particular, it follows that Voo 7^ 5q. Moreover, using ()3.2()( 
I3.21|) for z^oo and the equality in 1)3. 24() again, we also have 

-H''°°(0) > . 



On the other hand, — 11*^(0) is nothing but the expectation of v with respect to a positive, 
continuous function C{v) = 'Yld=2''^i\^i ~ u{v)\" . Hence 

- lim H'^J^O) := lim f C{v)vj{dv) > [ C{v)voo{dv) := -H^°°(0) . 
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This, together with (|3.23|) . imphes that 



eai^oo) > lim e^(i/j) . (3.25) 

Again, i^j is a maximizing sequence for hence there is an equahty in (|3.25() . so = eo-(i^oo)- 
By the second part of Lemma [3.31 we obtain the existence of € for which ea-{i^oo) = 

References 

[Am] Ambrosio, L, Lecture notes on Pptimal Transport problems, Lect. Notes in Math., 
"Mathematical aspects of evolving interfaces" (CIME seris Madeira (PT) 2000, 1812, 
P. Colh and J.F. Rodrigues, Eds., 1-52, 2003 

[AGK] J.E. Avron, O ; Gat and O. Kenneth, Optimal swimming at low Reynolds numbers, 
Phys. Rev. Lett., 93, # 18, 2004. 

[CHK] Chipot, M ; Hastings, S and Kinderlehrer, D: Transport in a molecular system, M2AN 
Math. Model. Num. Anal., 38 , #6, 1011-1034, 2004 

[CKK] Chipot, M; Kinderlehrer, D and Kowalczyk, M: A variational principle for molecular 
motors, dedicated to Piero Villagio on the occasion of his 60 birthday, Mechanica 38, 
#5, 505-518, 2003 

[DKK] Dolbeault, J; Kinderlehrer, D and Kowalczyk, M:Remarks about the flashing rachet, 
Partial differential equations and inverse problems, Contemp. Math. 362, 167-175, 
2004 

[E] Evans, L.C.: Weak convergence methods for no nonlinear partial differential equations, 
CBMS, 74, 1990 

[W] Wolansky, G: Rotation numbers for measure-valued circle maps, J. D'analyse, to appear. 



11 



